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Abstract 

We study the universal critical behaviour near weakly first-order phase transitions for a 
three-dimensional model of two coupled scalar fields - the cubic anisotropy model. Renorma¬ 
lization-group techniques are employed within the formalism of the effective average action. 
We calculate the universal form of the coarse-grained free energy and deduce the ratio of 
susceptibilities on either side of the phase transition. We compare our results with those 
obtained through Monte Carlo simulations and the e-expansion. 


A quantitative description of weakly first-order phase transitions (which are typically fluctu¬ 
ation driven) is necessary for a wide range of problems. Apart from their relevance to statistical 
systems (such as superconductors or anisotropic systems, like the one we consider below), weakly 
first-order phase transitions often appear in the study of high-temperature field theories. As a re¬ 
sult they have important implications for cosmology. An example is provided by the electroweak 
phase transition, which is weakly first order for Higgs boson masses between approximately 40 
and 80 GeV ||, |2|. The calculation of the baryon number that can be generated during the 
electroweak phase transition requires a precise determination of its properties. 

A simple example of a system with an arbitrarily weakly first-order phase transition is the 
cubic anisotropy model § 1- In field-theoretical language, this corresponds to a theory of 
two real scalar fields 4> a (a = 1,2) in three dimensions, invariant under the discrete symmetry 
(1 <-> — 1,2 <-> — 2,1 2 ). It can also be considered as the effective description of the four- 

dimensional high-temperature theory with the same symmetry, at energy scales smaller than 
the temperature i 10- It is relevant for the phase transitions in multi-Higgs extensions of the 
Standard Model. 

Recently, the properties of the weakly first-order phase transitions in this model have been 
studied in detail, within the e-expansion Q or through Monte Carlo simulations [|9]. In par¬ 
ticular, universal amplitudes have been computed, which describe the relative discontinuity of 
various physical quantities along the phase transition, in the limit when the transition becomes 
arbitrarily weakly first order. A discrepancy has been observed in the predictions for the uni¬ 
versal ratio of susceptibilities x+/X— on either side of the phase transition obtained through 


Monte Carlo simulations or the e-expansion [1C]. More specifically, the Monte Carlo simulations 
predict x+/x~ = 4.1(5), while the first three orders of the e-expansion give x+/x~ = 2-0,2.9, 2.3 
respectively. 

In this letter we present an alternative approach that employs the exact renormalization 


group [111 . It is based on the effective average action T /. [ 121 , which is a coarse-grained free 
energy with an infrared cutoff. More precisely, T k incorporates the effects of all fluctuations with 
momenta q 2 > k 2 , but not those with q 2 < k 2 . In the limit k —> 0, the effective average action 
becomes the standard effective action (the generating functional of the 1PI Green functions), 
while at a high momentum scale (of the order of the ultraviolet cutoff) k = A —> oo, it equals 
the classical or bare action. It is formulated in continuous (Euclidean) space and all symmetries 
of the model are preserved. The exact non-perturbative flow equation for the scale dependence 
of T fc takes the simple form of a renornralization-group-improved one-loop equation 111 
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The trace involves a momentum integration and summation over internal indices. The relevant 

( 2 ) 

infrared properties appear directly in the form of the exact inverse average propagator TL , 
which is the matrix of second functional derivatives with respect to the fields. There is always 
only one momentum integration - multiloops are not needed - which is, for a suitable cutoff 
function R k (q 2 ) (with R k ( 0) ~ k 2 , R k (q 2 —> oo) ~ e~ q ), both infrared and ultraviolet finite. 

The flow equation ( 0 ) is a functional differential equation, and an approximate solution 
requires a truncation. Our truncation is the lowest order in a systematic derivative expansion 
of T fc |4, “ 
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The wave-function renormalization is approximated by one ^-dependent parameter Z k . The 
truncation of the higher derivative terms in the action is expected to generate an uncertainty 
of the order of the anomalous dimension 77 . (For 77 = 0 the kinetic term in the P-dependent 
inverse propagator must be exactly proportional to q 2 both for q 2 —» 0 and q 2 —> 00 .) For the 
model we are considering, 77 ~ 0.035 and the induced error is small. A similar truncation has 
been employed for the calculation of the equation of state for the O (IV)-symmetric scalar theory 
fll 6 || . The result is indeed in agreement with those obtained through alternative methods with 
an accuracy of order 77 . 

The evolution equation for the potential results from the substitution of eq. ([ 2 ]) in the exact 
flow equation (|l ; ). The fixed-point structure of the theory is more easily identified if we use the 
dimensionless renormalized parameters 

Pa = Z k k~~ l p a Uk(Pl,P2)= k~ 3 U k (pi, P2)- (3) 


The evolution equation for the potential can now be written in the scale-independent form 

d 11 

■^u k (pi,p 2 ) = -3 u k + (1 + ri)(piui + p2U 2 ) - ^Ll(m\) - ^L 3 0 (ml), 


( 4 ) 


where t = \n(k/A). The anomalous dimension 77 is defined as dlnZ k /dt = —77 and can be 


computed starting from the exact flow equation 112, 14]. For the model we are considering and 
the part of the phase diagram we shall study, it is constant 77 ~ 0.035 to a good approximation 
[7j]. The quantities rh 2 2 are the eigenvalues of the rescaled mass matrix at the point (p \, p 2 ) 
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1,2 = x \ u i + u 2 + 2unpi + 2u 22 p2 ± (ui - u 2 + 2uupi - 21122 P 2 ) + 16mi 2 pi/5 2 
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and we have introduced the notation u\ = du k /dpi, u\ 2 = d 2 u k /dpidp 2 , etc. The function 
Lq(w), as well as the functions L\{w) = — dLl(w)/dw , L 3 n+1 (w) = — 1 /n dLf l (w)/dw for n > 1 
that we shall encounter in the following, are negative for all values of w. Also \L^(w)\ are 
monotonically decreasing for increasing w and introduce a threshold behaviour in the evolution. 
For large values of rh 2 the last two terms in eq. (||) vanish and the evolution of U k stops. The 
above functions have been extensively discussed in refs. [ 


0 - 


The initial condition for the integration is provided by the bare potential, which is identified 
with the effective average potential at a very high scale k = A. We use a bare potential of the 
form 

ua(pi,P 2 ) = {(pi - k a ) 2 + {p 2 ~ «a) 2 } + (1 + x\)X\pip 2 ( 6 ) 


and Z\ = 1. 

The phase structure of the theory has been discussed in detail in refs. 10 , 10 - The phase 
diagram has three fixed points, which govern the dynamics of second-order phase transitions. 
They are located on the critical surface separating the phase with symmetry breaking from 
the symmetric one. The most stable of them corresponds to a system with an increased 0(2) 
symmetry. It can be approached directly from a bare action given by eq. ([(j|) with x\ = 0. 
The other two are Wilson-Fisher fixed points, corresponding to two disconnected ^-symmetric 
theories. One of them can be approached from a bare action with x\ = — 1 , while the second 
requires x\ = 2 Q Flows that start with —1 < x\ <2 and near the critical surface eventually 

1 A redefinition of the fields demonstrates that this choice corresponds to two disconnected ^-symmetric 

theories Jt|. 
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lead to the 0(2)-symmetric fixed point. For x\ > 2 or x\ < —1 the evolution leads to a region 
of first order phase transitions. If x\ is chosen slightly larger than 2 or slightly smaller than —1 
the phase transitions are weakly first order. As we demonstrate in the following, the evolution 
first approaches one of the fixed points before a second minimum appears in the potential. 

The partial differential equation ( 0 ) with rj = 0.035 can be integrated numerically through 


a generalization to the case of two fields of the algorithms presented in ref. [17] for the case of 
one background field. This approach is straightforward but requires excessive computer power. 
An alternative solution relies on an approximation that simplifies the form of the potential. We 
are interested in the case where the minima of the potential are located on the two axes (this 
requires x\ > 0). We concentrate on the pi axis and consider the region p 2 -C 1. The potential 
can be expanded as 

Uk(puh) = Vk(pi) + fk(pi)P2 + ( 7 ) 

Substituting the above expression in eq. (|4|) and setting p 2 = 0 results in an evolution equation 
for Vk(pi)- Notice that terms with powers of p 2 higher than the first do not affect the form of 
this evolution equation. 

We approximate the function fk(pi) by a fourth-order polynomial in p\. This is expected 
to be a good approximation as the physical behaviour of interest (scale-invariant form of the 
potential, appearance of a minimum at the origin) is observed for p\ < 1 (see fig. 1). Near the 
origin we have 

fk(Pi) - ml + (1 + x s ) A s pi, (8) 

with ml = v\(p\ =0), A s = vu(pi = 0). The form of the first term in the right-hand side is 
imposed by the 1 2 symmetry. The evolution of the ^’-dependent term x s can be calculated 

from a truncated form of the evolution equation at p\ = p 2 = 0 [JTj 


(x s + l)x s (x s - 2)X s Ll(m 2 s ). 


dx s 1 

dt 8-7T 2 ' 

Near the minimum of the potential at pi = k ^ 0 we have 


fk(Pl) - x 0 X 0 K + (1 + xq) A 0 (pi - k) , 


(9) 


( 10 ) 


with Ao = vn(pi = k). Again the form of the first term in the right-hand side is imposed by the 
1 <-> 2 symmetry. The evolution of the /c-dependent term x$ can be calculated from a truncated 
form of the evolution equation at p\ = k, p 2 = 0 P] 


dx o 
dt 


1 6x 0 + ^ 
87T 2 K 1 — 


|7'i(2Aok) — L \ (xoAok)| + ^2* t o^o |9L2(2Aok) + (1 + xo) 2 T2(xoAok)| . 

( 11 ) 

Matching the two expansions of eqs. (| 8 |) and (1C) fixes the coefficients of the fourth-order 
polynomial. We shall check the validity of this approximation in the following. 

The numerical integration of the partial differential equation (|4]) is now feasible, because 
we have only two independent variables (t,pi) for the function vk (as well as the /c-dependent 
parameters x s , xq). The expected numerical accuracy of the integration is estimated to be of 
order 3% |L7]]. 

We are interested in the region x > 2. The region x < — 1 can be mapped on it through a 
redefinition of the fields @. In fig. 1 we present the evolution of v\ = dvk/dpi, starting at a 
very high scale k = A with a bare potential given by eq. (p). All dimensionful quantities are 
normalized with respect to A. The initial coupling Aa is chosen arbitrarily, while the minimum 
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aca of the bare potential is taken very close to the critical value K cr that separates the phase 
with symmetry breaking from the symmetric one. For |5ka| = |«a ~ H cr \ 1 the system 
spends a long “time” t of its evolution on the critical surface separating the two phases. The 
initial value x\ (which determines the initial conditions x s (A) = xo(A) = x\ for eqs. ©, ©) 
is taken slightly larger than the fixed-point value 2. After the initial evolution (dotted lines) 
the potential settles down near the Wilson-Fisher fixed point (solid lines). The fixed point 
of eqs. (D, (0) at x s = xq = 2 is repulsive and x s , xq eventually evolve towards larger 
values. This forces the potential to move away from its scale-independent form (dashed lines). 
At some point in the subsequent evolution the curvature of the potential at the origin v\{p\ = 0) 
becomes positive. This signals the appearance of a new minimum there, and the presence of a 
fluctuation-driven first-order phase transition. The evolution of the potential after it moves away 
from the fixed point is insensitive to the details of the bare potential. It is uniquely determined 
by |<5ka| = |^a — n C r | 1 and 5x a = xa — 2 <C 1, and, therefore, displays universal behaviour 


In the limit k —> 0 and in the convex regions, the rescaled potential uk grows and eventually 
diverges in such a way that Uk becomes asymptotically constant, equal to the effective potential 
U = Uq. This is apparent in fig. 2, where the potential along the p\ axis is plotted. The 
evolution of the non-convex part of the potential (between the two minima) is related to the 
issue of the convexity of the effective potential. This part should become flat for k —> 0 [18, |Tg| l . 
The approach to convexity is apparent in fig. 2, even though we have not followed the evolution 
all the way to k = 0. The nucleation rate during the decay of the unstable minimum can be 
calculated in terms of the potential at an appropriate coarse-graining scale. It is exponentially 
suppressed by the free energy of the dominant tunnelling configuration. The scale k must 
be chosen so that the pre-exponential factor, arising from fluctuations around the dominant 
configuration, is small. An appropriate scale kf for the definition of a universal coarse-grained 
potential is determined by [19] 


k 2 

K f 


(< d 2 U kf /d </>?) 




kj 


= 0 . 01 , 


( 12 ) 


where 4>\ ma x corresponds to the top of the barrier. The scale at which we have stopped the 
evolution in fig. 2 is near kf. 

The relative magnitude of \Sk\\ = \ka — H cr \ 1 and 5x a = sa - 2 < 1 results in different 

types of evolution. For the type of behaviour depicted in figs. 1 and 2 one must take |<5 /£a | <S Sx a- 
In the opposite limit, |<5 ka| Sx\, the system leaves the critical surface before x s , xq evolve 
away from the fixed-point value 2. As a result, a second minimum never appears at the origin. 
Instead, the only minimum of the effective potential U$ is located either at zero (symmetric 
phase) or away from it (phase with symmetry breaking), depending on the sign of 6k\. The 
resulting phase transition is second order and occurs for = 0. 

We are interested in the universal ratio of susceptibilities x+/x~ on either side of the 
phase transition. This ratio depends on the value of \5k&\/5x\. For \5k\\/5x\ > 1 we have 
x s (k),xo(k) — 2 during the whole evolution. The universal quantities characterizing the second- 
order phase transition are determined by the Wilson-Fisher fixed point. We calculate x+/x~ 
by integrating the evolution equation and evaluating cPUo/dcpi = x 1 the minimum, for 

2 Despite the presence of 1 — xo/2 in the denominator in the right-hand side of eq. (|i|), this equation has a 
fixed point at xo = 2. 
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5k a = Te with e <C 1. We obtain x+/x~ = 4.1. This value should be compared with the 
value x+/x~ = 4-3, calculated with the same method for the one-field, ^-symmetric theory 
(the Ising model) without any polynomical approximations for the form of the potential [p^|. 
The difference is due to the approximation of the function fk(pi) in eq. by a fourth-order 
polynomial. An error of this magnitude is typically induced by the truncation of the potential 

]. We have verified this conclusion by approximating the function 
without using eq. ©• This cruder approximation results in \+/X- = 4-0. 


to a polynomial form [14 
fk(pi) by eq. (f 


The evaluation of the same quantity through the e-expansion or an expansion at fixed di¬ 
mension gives x+/x~ = 4.8(3) pO], whereas experimental information gives x+/X— = 4.3(3) 


[20]. The difference with the results of our method is due to the omission of higher derivative 
terms in the effective average action of eq. (jp). This confirms our expectation, discussed earlier, 
that the error induced by this approximation is determined by the anomalous dimension, which 
is small for our model (77 ~ 0.035). 

For \8k\\/5x\ -C 1 the potential develops a second minimum at the origin during the later 
stages of the evolution. The phase transition is approached by fine tuning 5k\, so that the 
two minima have equal depth for k = 0. In fig. 3 we plot the difference in energy density A U 
between the minimum away from the origin and the one at the origin, as a function of the scale 
k, for four values of 5k\. The numerical integration of the evolution equation is difficult for 
k —► 0 because of the singularity structure of L 3 n (w) [[T^, 19|. As a result the curves of fig. 3 
must be extrapolated to k = 0. Line (a) corresponds to a system in the symmetric phase, line 
(d) to one in the phase with symmetry breaking. Lines (b) and (c) correspond to a system very 
close to the phase transition. The evolution of x+/X- f° r the same parameters is depicted in 
fig- 4. 

The extrapolated value of x+/X- for k = 0 near the phase transition is expected to lie in 
the interval (1.5, 2). We have used numerical fits of various curves for A U and x+/x~ i n the 
vicinity of the phase transition in order to perform the extrapolation to k = 0. The expected 
value for x+/x~ is 1-7. 

The main sources of error in our calculation are the derivative expansion of eq. (^) and the 
approximation of fk(pi) in eq. (j7|) by a fourth-order polynomial in fj\. We argued earlier that the 
error induced by the derivative expansion is related to the small anomalous dimension 77 ~ 0.035. 
We also checked above the error induced by the polynomial appoximation of the potential. For 
the second-order phase transition, two polynomial approximations give x+/X- = 4.0 and 4.1, 
while no polynomial approximation gives x+/X- = 4-3 in the order of the derivative expansion 
we are working [[T(|. The numerical integration and the extrapolation used in the calculation 
of x+/x - induce smaller errors. The e-expansion or an expansion at fixed dimension gives 
X+/x~ = 4.8(3) [^0|, whereas experimental information gives x+/x~ = 4.3(3) |2(J. In order to 
be conservative, we use as a maximum total error the difference between our result x+/X— = 4.1 
for the case of the second-order phase transition and the central value x+/x~ = 4-8 of the result 
from the e-expansion. 

Our result 

X+/X- = 1-7(7) (13) 


is depicted by a horizontal band in fig. 4. It is in good agreement with the predictions of the 
e-expansion (x+/x~ = 2.0,2.9, 2.3 for the first three orders) HU- It disagrees, however, with 
the lattice result presented in refs. [||. (x+/x~ = 4-1(5)). The values favoured by the lattice 

calculation correspond to a theory well into the phase with symmetry breaking in our calculation 
(line (d) in figs. 3 and 4). However, it should be noted that, according to the authors of ref. 
fllOll , the error of the lattice result “should be taken with a grain of salt”. 
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In summary, we presented a method that can provide a quantitative description of the 
universal behaviour near weakly first-order phase transitions. It is based on the calculation of 
a coarse-grained free energy through an exact flow equation. Fixed points in the evolution, 
the appearance of new minima in the potential, and the universal properties of the resulting 
fluctuation-driven first-order phase transitions can be studied in detail. We calculated the 
universal ratio of susceptibilities x+/x- on either side of the first-order phase transition in 
a two-scalar model . 
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Fig. 2: The potential along the p± axis, as the coarse-graining scale is lowered. 
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Fig. 3: The difference in energy density between the two minima of the potential, as a function 
of the coarse-graining scale, near the phase transition. Line (a) corresponds to a system in 
the symmetric phase, line (d) to one in the phase with symmetry breaking. Lines (b) and (c) 
correspond to a system very close to the phase transition. 
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